#----------------------------------
#----------------------------------
# WALAIR++ Input File
# Alastair Chester
# PhD Paper 2
# 100kg TNT Charge at 25m on gauge block
# 2_7_6 - 100mm
#----------------------------------
#----------------------------------


# GPU Command - Only needed on multiple GPU connected systems
# Device 0


#----------------------------------
#Chuncks - Partition domain if model too big for GPU memory - default is 1
#CHUNKS 1


#----------------------------------
#Simulation end time in seconds

TIME 0.1

#----------------------------------
#CFL1D - recommended 0.5 unless JWL is used then potentially lower 

CFL1D 0.5


#----------------------------------
#CFL2D - recommended 0.5 unless JWL is used then potentially lower 

CFL2D 0.5


#----------------------------------
#CFL3D - recommended 0.5 unless JWL is used then potentially lower 

CFL3D 0.5


#----------------------------------
#Time step control - optional command, generated internally but can be overridden here
#TSTEPCONTROL abs 1.0e-8 


#----------------------------------
# Run Type - choose the combination of domains and remaps 11, 22, 33, 12, 13, 23 or 123 for 1D-2D-3D remap

RUNTYPE 123


#--------------Charge and Air Definitions-----------------
#Set ambient conditions if differnt to defaults pressure = 101325 Pa, density = 1.225 kg/m3 and γ = 1.4.

#AMBIENT 1.09e5 1.4 1.1


#----------------------------------
#Set Charge Mass in kg

CHARGEMASS 100


#----------------------------------
#Set Charge Properties density energy
#TNT values from manual

CHARGEPROP 1600.0 4.52e6 


#----------------------------------
#Set the 1D Domain size, domain size in m

DOMAIN1D 1.2


#----------------------------------
#Set the 1D Cell number - Note this is the number of cells across the charge radius not the domain.

CELLS1D 50


#----------------------------------
#Set the 2D Domain size, r Zlower Zupper

DOMAIN2D 24.9 0 26.1


#----------------------------------
#Set the 2D Cell number - domain radius number, vertical cell size is set to same as radial.

#so for 30mm cell size on 24.9m 2D domain

CELLS2D 830


#----------------------------------
#Set 2D Boundary Conditions brl bru bzl bzu, 1 transmissive, -1 reflective. - If used 2D model will run to endtime. bzl is reflective by default and other boundaries and transmissive and terminate.

#BOUNDARY2D 


#----------------------------------
#Set the 3D Domain size xmin ymin zmin xmax ymax zmax 

DOMAIN3D -25 -10.75 0 5.45 10.75 9.95


#----------------------------------
#Set the 3D Cell numbers - number of cells for each dimension nx ny nz.

#so for 100mm cell size on above 3D domain

CELLS3D 305 215 100


#----------------------------------
#Set 3D Boundary Conditions xl yl zl xu yu zu, 1 transmissive, -1 reflective. - All boundaries transmissive by default

BOUNDARY3D 1 1 -1 1 1 1


#----------------------------------
#Auto domain sizing for 1D and 2D domains, based on 3D domain size, off by default (0)

#AUTODOMAIN 0


#----------------------------------
#Set Charge Location 3D coords x y z

CHARGELOC -25 0 1.2


#----------------------------------
#Optional - change z height from 2D to 3D domains

#REMAPZ 5


#----------------------------------
#If non-spherical charge sizes are required then can be entered here - see manual pg 57

#----------------------------------
#Secondary sympatheticly detonated charges can be included here see manual pg 64

#----------------------------------
#JWL can be included here - see manual pg. 65

#----------------------------------
#Afterburn see manual pg. 67

#----------------------------------
#setdet point for 2D and 3D charge definitions


#----------------------------------
#Checkpoint  - determine restart file write time frequency in seconds - optional

#CHECKPOINT 0.1


#----------------------------------
#Restart checkpoint for 1D in m, require command for each file

#RCHECKPOINT1D 0.1
#RCHECKPOINT1D 0.2


#----------------------------------
#Restart checkpoint for 2D in m, require command for each file

#RCHECKPOINT2D 0.1
#RCHECKPOINT2D 0.2


#----------------------------------
#Restart checkpoint for 3D in m, require command for each file

#RCHECKPOINT3D 0.1
#RCHECKPOINT3D 0.2


#----------------------------------
#Call a Restart checkpoint for 1D e.g.

#RESTART1D 1D_walairAnalysis_0.4.rstrt


#----------------------------------
#Call a Restart checkpoint for 2D e.g.

#RESTART2D 2D_walairAnalysis_0.4.rstrt


#----------------------------------
#Call a Restart checkpoint for 3D e.g.

#RESTART3D 3D_walairAnalysis_0.4.rstrt


#----------------------------------OUTPUT----------------------------------
#Set VTK slice or volume output frequencies - title slice_freq surf_freq vol_freq (-ve to turn off)

OUTPUT 2_7_6


#----------------------------------
#Output control - choose what should be output e.g. pressure density etc.

OUTPUTCONTROL pressure 1 overpressure 1 max_pressure 1 max_overpressure 1


#----------------------------------
#Human Injury factors can be included here - see manual


#----------------------------------
#Choose VRTK slice(s) x location

#VTKSLICEX 0.6
#VTKSLICEX 0.5
#VTKSLICEY
#VTKSLICEX


#----------------------------------
#VTK Box output xL yL zL xU yU zU - Optional - if larger than domain will revert to domain dimensions

VTKBOX -25 -13.75 0 13.45 13.75 13.95


#----------------------------------
#Are user obsticles required in output VTK files default (1) off

USERSURFOUTPUT 0 


#----------------------------------
#Pressure gauges in 2D

#PRESSUREGAUGE2D 2.1 1.3


#----------------------------------
#Pressure gauges in 3D




#----------------------------------
#Pressure file definition - read in pressure gauges from external file, data output file given for each input file

PRESSUREFILE ./pGauges100mm.txt


#----------------------------------
#Set pressure write frequency in computational steps, default is 1

#PRESSUREFREQ 10

#----------------------------------
#Ray tracing gagues can be included here see manual pg 74


#----------------------------------DEFINE Obsticles------------------------------
#Define simple obsticle shapes using OBSRECT3D, OBSCYL3D, OBSSPHR3D, OBSNEGRECT3D, OBSNEGCYL3D and OBSNEGSPHR3D

#OBSRECT xmin ymin zmin xmax ymax zmax

OBSRECT3D    0 -3.75 0 3.45 3.75 3.95
OBSNEGRECT3D 0.225 -3.75 0.225 3.225 -1.75 3.725
OBSNEGRECT3D 0.225 1.75 0.225 3.225 3.75 3.725
OBSNEGRECT3D 0.75 -1.75 0 3.45 1.75 3.95



#----------------------------------
#Can also load in stl files here
